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ABSTRACT 


The objective of this thesis is to find out the 
parametric instability regions and steady state response 
characteristics in the neighbourhood of main parametric 
resonance of a vertical bean which has a concentrated mass 
and is restrained at the upper end and simultaneously subje- 
cted to a harmonic axial displacement excitation at the 
bottom end. This system is a dynamic model of an engine- 
valve mechanism where the displacement excitation is produ- 
ced by a cam. The beam under this loading is considered to 
be executing lateral vibration with large amplitude which 
results in longitudinal displacement of any section and the 
axial displacement of the end load. The equation of motion 
for this lateral vibration of the beam is deduced. This has 
non-linear terms because of large curvature and longitudinal 
inertia of the beam elements and the concentrated mass. Two 
mode approximation for the lateral vibration of the beam is 
made and Galerkin's method is applied to reduce the governing 
partial differential equation to two coupled non-linear 
ordinary differential equations with periodic coefficients. 
These equations are linearised and then the instability 
regions are found out by the use of Hills determinants. The 
principal parametric instability region of these equations 
is compared with the corresponding region obtained from 



the linearised equation of one mode approximation. 

For the steady state response characteristics 
of one mode approximation in the neighbourhood of its 
principal parametric resonance the methods of (l) harmonic 
balance (2) slowly varying parameters and (3) multiple 
scales have been used and these results are mutually compa- 
red. The effects of various parameters like the ratio of 
longitudinal spring stiffness to the transverse tip stiffness 
of the beam,, ratio of the concentrated mass to the beam mass 
and the non-dimensional amplitude of the displacement excita- 
tion ( with respect to beam length) on the response curve 
are investigated. 
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CHAPTER 1 

INTRODUCTION 


In contrast with the forced oscillation of a 
single or multi-degree of freedom system (mechanical or 
electrical) in which excitation appears as inhomogenei- 
ties in the governing differential equation, in the 
parametically excited systems the excitation appears as 
the coefficients in the governing differential equation. 

So for these systems one is led to differential equations 
with time varying coefficients. Since the excitations 
when they are time independent appear as parameters In the 
governing differential equation, this type of excitation 
is called parameteric excitation. Moreover, in constrast 
with the external excitation in which a small excitation 
cannot produce large response unless the frequency of exci- 
tation is close to one of the natural frequencies of the 
system, a small parameteric excitation can produce a large 
response when the frequencies of the excitation are other 
than the natural frequencies of the system. Moreover, 
parameteric instability occurs over a region of a parameter 
space consisting of excitation amplitude and frequency and 
not at discrete points. This instability can correspond 
to large displacement in the direction normal to the longi- 
tudinal excitation of finite amplitude. 
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Parametrically excited systems, reducible to the 
standard Mathieu Equation having two parameters (6 and E ) , 
have some boundary curves in the parameter plane (d,£) which 
separate the stable and unstable regions. The stable and 
unstable regions are of interest because one wants to avoid 
the unstable regions so that there is no possibility of very 
large response of the system. 

There are a number of analytic techniques availa- 
ble for the determination of the boundaries separating 
stability from instabilities. These techniques can be 
divided broadly into following three classes- 

(i) Hills Method of Infinite Determinant 

(ii) Perturbation Methods 

(iii) Liapunov’s Theory. 

The Hill’s method of infinite determinants is 
extensively used for a single degree of freedom system and 
has been recently applied to multidegree of freedom systems 
by a number of researchers [1,2]. The second class consists 
of perturbation methods [16] which assume that tho variable 
coefficient terms are small in some sense. It is also 
necessary to assume a series solution of the differential 
equation in powers of a small dimensionless parameter G and 
make sure that it is a uniformly valid exponsion which 
means that the higher order terms in S should provide a 
small correction to the lower order terms of S . 
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The third class uses Liapunov's theory [15]. This approach 
is limited by the ability to find out a suitable Liapunov 
function. 

The influence of damping on the instability boun- 
darios was discussed by many researchers [3,4,5]. In most 
cases, damping force has stabilizing effect but Valeev [6] 
showed that for certain combination resonances the damping 
force may alter a stable state into an unstable one. 

Moreover, Stevens [8] showed that some viscoelastic materials 
are destabilizing. 

By deleting all the non-linear terms in the gover- 
ning equation of the system, one arrives at a differential 
equation, called the corresponding linear system. It plays 
a key role in the analysis of weakly non-linear systems. One 
can write the solution of a non-linear system in the form of 
a perturbation series which may bo considered as a develop- 
ment in the neighbourhood of the solution of the corresponding 
linearised equation. There are a number of ways in which 
this perturbation can be effected. One starts with straigh- 
tforward expansion, which is not cniformaly valid. So 
several modifications of the straight forward procedure are 
to be done so as to lead to uniformly valid expansion. One 
modification of it is the Lindsted-Poincare Method where a 
new time scale is introduced to permit the frequency and 
amplitude to interact and there are provisions to avoid 
secular terms ( which give, unbounded solution). But this 
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method can not treat damped systems conveniently. So 
method of multiple scales is used as a further modification 
of Lindsted -Poincare's Method, The underlying idea of 
the method of multiple scales is to consider the expansion 
respresenting the response to a function of multiple indepen- 
dent variables or scales, instead of a single variable. The 
steady state vibration can also be analysed by the method of 
hormonic balance and by method of slowly varying parameters. 

BRIEF REV I £•••; OF AVAILABLE LITERATURE : 

Parametric instability was first observed by 
Michael Faraday in 1831 [7], when he found that surface 
waves in 'a fluid-filled cylinder under vertical excitation 
exhibited twice the period of excitation itself. Melde [7] 
in 1859 found that though a longitudinal force is applied 
on a string tied between a rigid support and a heavy tuning 
fork, it exhibits a lateral response. Beliaev in 1924 first 
considered the parameteric instability of a column pinned 
at both ends and subjected to a time dependent axial force 
[7;19]. He reduced the governing equation of motion to a 
standard Mathieu-Hill equation and showed the stability and 
instability regions in the parameter space. Later many other 
scientists extended his work, Bolotin [3] considered various 
kinds of structural elements and discussed their dynamic 
stability. He derived a modified equation of motion for a 
column considering the initial curvature and longitudinal 
inertia of the column elements and the attached end mass. 
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Evensen and Evan- Iwanowski [9] and Handoo and Sundararajan [10] 
investigated the parametric instability regions of a perfect 
column with an attached end mass both analytically and experi- 
mentally. They considered the longitudinal inertia of the 
beam and the attached mass. Carlson [11] considered a tensi- 
oned bar with initial curvature and an elastic spring restra- 
int at an end which included a lumped mass. He found the 
steady state response characteristics by the method of har- 
monic balance. Sato and Saito [12] investigated the steady 
state response of parametric vibration of a simply supported 
horizontal beam, carrying a concentrated mass which is not 
at the end under the influence of gravity and periodic 
longitudinal exciting force. They included non-linear terms 
arising from moderately large curvatures, longitudinal iner- 
tia of the beam and concentrated mass and rotatory inertia of 
the concentrated mass in the equation of motion. They appl- 
ied harmonic balance method to investigate the dynamic 
response. Saito and Koizumi [13] investigated the steady 
state response of parametric vibration of a simply supported 
horizontal beam with a concentrated mass at one end and 
subjected to a peroidic axial displacement excitation at the 
other end under the influence of gravity. They considered the 
non-linear terms arising from longitudinal inertia of concen- 
trated end mass and beam clement. They also used harmonic 
balance method for finding the steady state behaviour. 

Gurgoze [14] found out the principal parametric instability 
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region and steady state response of a simply supported 
vertical beam by Mettler's method and method of harmonic 
balance respectively. The beam carries a concentrated end 
mass and is restrained at one end and subjected to a perio- 
die axial displacement excitation at the other end. He 
considered one mode approximation of its lateral vibration 
and included non-linear terms arising from moderately 
large curvature and longitudinal inertia of the concentrated 
mass and beam element. 

SCOPE OF THE PRESENT WORK : 

In the present work, a simply supported vertical 
beam whose upper end is restrained and carries a concentra- 
ted mass and the lower end is subjected to a harmonic exci- 
tation force is considered. The large curvature deflection 
theory gives a partial differential equation for its lateral 
vibration. This equation has non-linear terms due to large 
curvature and longitudinal inertia of the concentrated mass 
and the beam elements. Here two mode approximation of the 
lateral vibration of the beam is considered and the gover- 
ning partial differential equation of motion is reduced 
to two non-linear ordinary differential equations. After 
linearising these equations the instability boundaries are 
found out making use of Hills determinants. Taking one 
mode approximation of lateral vibration of the beam the 
instability boundaries and the effects of damping, mass 
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ratio and stiffness ratio on it are found out and hence, 
the result obtained by Gurgoze is ’reproduced. Taking one 
mode approximation the backbone curves are found out by 
method of harmonic balance, method of slowly varying 
parameters and method of multiple scales. 

Chapter 2 deals with the formulation of the 
system i.e. to derive the equation of motion of the system. 
Chapter 3 deals with the instability boundaries which are 
found out by the method of Hills determinants . These 
boundaries are found out for both ono mode approximation 
and two mode approximation. Chapters 4-6 deal with the 
method of hormonic balance, the method of multiple scales and 
method of solwly varying phase and amplitude to find out 
the steady state response. The data of reference [14] is 
used to get the backbone curves and these curves are plo- 
tted. The effects of mass ratio, stiffness ratio and 
excitation amplitudes (non-dimensional ) on the response 
curve are found out. And the results from all these three 
methods are compared with each other. In fact the results 
from method of harmonic balance are reproduced as it was 
done by Gurgoze [14], 
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CHAPTER 2 

FORMULATION 

A simply supported vertical beam of uniform 
cross-section, carrying a concentrated mass at one end and 
subjected to a periodic axial displacement excitation at 
the other end is considered, Fig. 1. The effects of shear 
deformation and rotary inertia of the beam are neglected 
as beam-thickness is small in comparison to the length. 
Bernoulli-Eulor theory is used for large lateral deflection 
of the beam and as such the coordinate system chosen is 
such that the distance of any section from the lower hinge 
is considered along the arc of the elastic curve and the 
lower hinge point is the origin. To construct the Bernoulli- 
Eulor equation, the moment of the spring force, longitudinal 
inertia forces of the concentrated mass and beam elements 
and lateral inertia force of the beam element are considered 
in the expression for bending moment at a cross section. 

These moments are calculated separately and are put in the 
Euler equation; differentiating the resulting equation 
twice with respect to the arc coordinate one gets the equa- 
tion of motion which is a non-linear partial differential 
equation (P.D.E.). Then one mode and two mode approximation 
for its lateral vibration are considered and using Galerkin's 
method, the governing P.D.E. is converted into non-linear 
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ordinary differential equation(s) (O.D.E.) which is/ are 
subsequently normalised to the standard form of Mathieu 
equation (s) . 

2 , 1 DESCRIPTION OF THE SYSTEM; 

In an internal combustion engine, the inlet and 
exhaust valves are operated by a cam system through engine- 
valve mechanism. The can lifts a cam follower connected to -f— 
push rod. The push rod actuates a lever, called a rocker 
arm, which depresses the valve stem and opens the valve. A 
typical valve actuating system for an air-cooled engine is 
shown in Fig. 2. 

This system is modelled, Fig. 1, such that the 
beam corresponds to push rod and mass M to the reduced mass 
of the rocker arm, valve and valve spring. The spring 
represents the equivalent stiffness of the valve spring 
and g (t) corresponds to the harmonic motion of the lower 
end of the push rod which is driven by a cam. The following 
assumptions are made: 

(i) The thickness of the beam is small in compari- 
son with the length of the beam. Therefore, the effect s of 
shear deformation and rotatory inertia are neglected. 

(ii) The extensional deformation of the beam is 


negligibl e 
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2.2 EQUATION OF MOTION ; 


From the Bernoulli-Euler theory, the bending 
moment at any section s is given, in the coordinate 
system shown in Fig. 1, by 


M. = El v 
o < 


ss 


(l~V 


1 

)~2 


( 2 . 1 ) 


where El is the bonding rigidity, v is the lateral displace- 
ment and subscript comma and S denotes differentiation with 
respect 1 to tho arc length S along the clastic curve. 

For tho parametrically oxcited beam shown in Fig.l, the 
bending moment at any section s can be written as 

M b = - F ( t ) v - -hM 2 -M 3 (2.2) 

where F(t) = sum of the spring and inertia forces excerted 
on the point mass. 

Mj- = bending moment at section s due to lateral 
inertia force 

M 2 = bending moment at section s due to longitudinal 
inertia force 

and = bending moment at section 5 due to damping 

force. 

The spring force at the upper end of the beam 
consists of preload on the spring F Q and the spring force 
due to the displacement of tho spring which is caused by 
ground excitation. To calculate the inertia force at the 
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concentrated mass, the displacement of the concentrated 
mass for any instant t is found out and it is differenti- 
ated twice with respect to time to get its acceleration. 
Therefore, the force F(t) is given by 


F(t) = F q +[g ( t)- (L-/ L Cos 0 d p) ] K f 

o 1 

L 9 

+ [g f tt~ / ( Cos +sin ^ * 0 tt^ u 

(2.3) 


where F Q 



L 

P 

and 0 


preload on the spring, 
spring stiffness and its coefficient 
is the displacement of upper hinge point 
of the beam, 

concentrated mass and its coefficient is 
the acceleration of upper hinge of the 
beam, 

length of the beam, 

arc length along the elastic curve 

angle between vertical axis and the tangent 

to the elastic curve. 


To find out the bending moment at any section S due 
to the lateral inertia force, one has to consider (3 along 
the beam which ranges from s to L. Then in this range one 
considers an element and finds out the lateral inertia of 
that element. To find out the total moment integrate the 
moment of the lateral inertia force over the entire range 
of (3 , Therefore 
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L (3 L 

M x = / \i v tt ( / Cos 0 drj) dp -R/ Cos 0 dp .... (2.4) 

s ' s s 

whore p = beam mass per unit length 

and Ft = reaction force at the upper end of the beam.Seo 
equation (2. 7 ) . 

Similarly, to find out the bending moment at section s due 
to longitudinal inertia force one considers p along the beam 
which ranges from s to L and is given by 

L p p 

M 9 = / pig ++ +(p~ / Cos 0 dr) ) .,] ( / Sin 0 dr)) dp 
z s o s 

... (2.5) 

where 

r) = arc coordinate along the elastic curve which 

ranges from S to p 

P 

and f Sin 0 drj = distance of the element from the vertical 
s 

axis. 

Similarly, the bending moments M ^ due to damping 
force distributed over s to L is given by 
L p 

Mo = / c v / Cos 0 dr)) dp ... (2,6) 

s ' s 

where & = viscous damping coefficient, 

r) = arc coordinate considered along the elastic 

curve which ranges from s to p 
p 

and / Cos 0 dr] = axial distance of the element from the 
s 

section s considered. 

To find out the horizontal reaction force R at 
the top end of the beam, take moment about the lower hinge 
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point of all the inertia forces and the damping force so 
as to obtain 

L 

R = [M ll(0) - M 2 (0) + M 3 (0)]/ / Cos 0 d p (2.7) 

where the first term in the expression for has been 
denoted as M-^. 

Now it is sufficient to express the terms 0, 0 0 ^ , 

Sin 0 and Cos 0 by their power series approximations because 
only upto third order non-linear terms are to be retained in 
the governing differential equation and hence 


0 


v 


P 


0 ,t= \pt 
0 , tt _ v , ptt 


Sin 0= 0=v 




Cos 0= 1— ^“ 2 “" 


1- i- V 2 

2 \{3 


Substituting these expressions into equations (2.3), (2.4) 
(2.5), (2.6) and (2.7), one obtains 

L l 2 

F(t) = F 0 + [g(t)-L+ / ( 1- ) d P 3 K f 

L 9 L 

+ [ g ftt - / dp - / v R v R++ dp] M ...(2,8) 


pt 


,P \ptt 


M, 


dll) dp 


R / ( 1 


1 2 \ 

- v fl ) dp 


. . (2.9) 
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M, 


/ 4 [g, tt + / (v >tJ t + V ,4 v ,ntt } dTl] 


• ( / 


^ v dr)) dp 


L p , 

M 3 = / C V t [/ ( 1- ~ V ,T)) d A d f 

S S 


.... ( 2 . 10 ) 
( 2 . 11 ) 


R 


r~ (A 4V + t [ / P ( 1“ b; V 2 ) dll ] dp 

j-p o » ^ L 0 ■<- tn 


L P P 

/ 4 [g ++ + f (v z „ + + v „ + + ) dii ] 

o ' o 


,T)t , r) ,r)tt y 
L 

’ J 
'P 


‘ ( drj) d 0) + T ( / cv , [/ (1- 

° > n d o * o 


'"^' 2 n ) dr)] dp) 
<*- ? ■ ) 


where L p = / (1- ™ v 2 ^ ) dp 


( 2 . 12 ) 

(2.13) 


projected length of the elastic curve on the 
vertical axis. 


Substituting the expressions (2,8) - (2.12) in 
equation (2.2) subject the equation (.2.13) and subsequently 
in the equation (2.1) and differentiating the resulting 
equation twice with respect to s one gets the following 
partial differential equation which contains all the principal 
non-linear terms up to the third order [ See Appendix -I for 
details ] , 


El 

r 1 2 

L , ssss 2 , ssss , s 

+ 3 

V „ 

V 


t s 

f s s 

+ 

\s \ss { cv ,t d f> + 

c \t 

U~ I v 2 s 


v + v° ] 
ss , sss , ss J 


ss 
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+ v v 
,s , ss 


- V 


„ V 

s . ss 



C g 


^ 2. \ 

S pv .j . ^ dp + p.v (l " 2 ^ ; s 

L P 2 

R+ ^ v ,ss l f 9,tt + { ( V„t + v ,n Vtt )dll] dp 

s o 

tt + { ( v fnt +v ,,i v ,ntt ) = 0 ••••(2.14) 


The boundary conditions for the simply supported 

beam are 

v= v = 0 at s = o and s =L ,...(2.15) 

f s s 

A solution to the differential equation (2.14) 
could bo sought in the form similar to that of lateral vibra- 
tion of beam with small amplitude which is given by 

oo 

v(s,t) = £ B (t) Sin SP. ( 2.16) 

n=l n L 

Experimental evidence shows that the first spatial 
mode is dominant. But here upto second mode is considered 
so as to investigate the effect of higher mode on the 
instability boundaries of the beam. 

2 . 3 O NE MODE APPROXIMATION ; 

One mode approximation for the equation (2.14) is 

v(s,t) =B(t) Sin ... (2.17) 

where B(t) represents a time dependent coefficient to bo 
determined from Galerkin's method. In this method, one puts 
the approximate solution in the loft hand side of the 
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equation (2.14) and the error function is orthogonalised 
with respect to the trial function . For one mode approxi- 
mation the error function is orthogonalised with respect 
to Sin only to obtain the following 0. D. E, which has 

been transformed to nondimensional eouation. 


zz + jaw z + ^ ^2 ^ OJ 01 “ K f 2 ® " M 22 a ,zz^ 

2 

+ -B < V ) “o + K fD 5 3 + M 46 B bA 

CO 7 


o 

3 -re c 

8"" pco 



B 


.. . (2.18) 


where, co is the angular frequency of the proscribed motion, 


z = 

cot, B = 
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“ L 
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2 0 
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K f2 

■re 2 K f 

_ 

' » K f4 

re 4 K f 

4pL ? M 22 ~ 

_2 r M 

71 LjlL + 

§]> 

M 46 

2 

re p 

2 *■ 

_2;M\ 

% ^pL^ 

aL f 3 n 

3 f 8 J 



and 

M 44 = 

re 2 r 2 
~2~ L ft 

M re 2 _ 3-i 

pL “ 3 8 J 


(2.19) 


Here, w = fundamental bending eigen frequency 
of the beam 

and co Q -p = fundamental bending eigen frequency of the same 
beam when compressed by the axial force F , 

Depending upon the cam profile there are different 
cam-displacement curves which arc periodic in nature. Also 
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in cam-f ollower mechanisms there can be loss of contact 
between the cam and the follower which depend upon its 
angular velocity and the com follower geometry. However, 
here it is assumed that the displacement excitation is a 
simple harmonic motion and there is no such loss of contact 
during motion. As such, the displacement is given by 

g(z) = B ( 1- Cos z) .... (2.20) 


Substituting (2.20) in equation (2.18), one gets 


B + 
, zz 


c 


B 


, z + ^ (i 2 ^ w ol “ 6 K f2^ + G ^ w 2 ~ M 22^ Cos 2 


+ 


w 

3 % : 

8 


[ (f) 


to 


+ K f4 ] B J + M a „ B B 


S2 


46 


iE“ ® 2 ' 5 , 2 + M 4 4 b 2 B, zz =0 


( 2 . 21 ) 


Equation (2.21) has seven terms and starting from the 
first these represent inertia, linear damping, linear stiff- 
ness and parametric excitation, the effect of non-linear 
curvature on the non- linear elasticity ( geometric non- 
linerity only), the non-linear longitudinal inertial force of 
tho beam itself, the non-linear character of damping and the 
non-linear inertia force of the concentrated mass, respectively. 


2.4 TWO MODE APPROXIMATION : 


For two mode approximation of tho response of the 
system, express v(s,t) as 
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v(s,t) = 


Bj_ ( t ) Sin 


+ B 2 (t) Sin 


2-rcs 

L 


... ( 2 . 22 ) 


where, the unknown functions B^(t) and S 2 (t) are evaluated 
by using Galerkin's method as discussed earlier. To this 
end B^,B 2 and g(t) are normalised with respect to L and the 
following two equations are obtained 


|iL 2 

Bj. + 

Clj2 B 
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2L 2 2 
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B2 



+ 'hi 

(%, 

• » » 

Bl , Bj; ) + 

' ^12 

1*2’%' 

* * 
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+ *13 
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b 2 , V 
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« • • • 

B 1 » B 2 ^ 
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... (2.23) 

and 
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f L + g 5 tt 
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K 




+ ^23 ^ B 1 * B 2 f B i * 

B 2 , b x , b 2 ) - 0 

... (2.24) 

where 

%(t) = 

B 1 (t) 

T“~ 7 


B 2 (t) = 

B 2 (t) 

L 


g(t) = 




where 
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and tho functions {jj^ are given by 
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^22 <M 2 'KP = +2Tc 4 K f L ) I 2 

Li 


+ [ 4it 4 ML + (4 +§2l) ^ 2 l 2 ] sr 2 fei 


+ B 2 B 2 [ 4h 2 ML + ( ” - ~— 


2 u 2 
) \xn 2 L 2 ], 
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4 H 2 °2 
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Now linearising the equations (2.23) and (2.24) and 
introducing the new time scale z = tot the following equations 
are obtained. 

* \,zz + 6 W B l,z + ( k - 2> B 1 

4 "ol 
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and 

q 2 n + e C* q EL + (1 + 
2,zz 



6K 


f 2 


to 


ul 



+ B ( ~ 4q 2 M22) Cos z B 0 - 

w ol 


40 e 
9 


q Cos z 8| = O 


( 2 . 25 ) 


where 


c* - 

and 

q = 


2|iBu) 0 i 


to 



These are two coupled linear equations with periodic 
coefficients, refered to as coupled iviathieu equations. 
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Fig. I Restrained beam with end mass 




Fig.2 Overhead -valve, rocker arm and push-rod 
assembly. 
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CHAPTER 3 

INSTABILITY ANALYSIS 

The present chapter deals with the stable and 

unstable regions in the E, q plane for two coupled Mathieu 

equations for two mode approximation of lateral vibration 

and 

of the beam* The parameters E^q refer to non-dimensional 
excitation amplitude and the ratio of excitation frequency 
to the bending eigen frequency of the beam respectively. 

The transition curves E = E (q, M*., S^) separating the 
stable and unstable regions for different parameters like 
mass ratio j\q. , stiffness ratio and damping S ^ are of 
prime interest. The Mathieu equations have 2 % - periodic 
coefficients, therefore, from Floquet theory it is concluded 
that for such values of E and q these Mathieu equations 
possess a periodic solution having either period 2 % or period 
4tc . Since these periodic solutions Bj(t) and B 2 (t) are regu- 
lar for all values of time it follows that these can be expre- 
ssed by appropriate Fourier series of 2 % and 4 % periods. 
Substituting these in the DE of (t) and Bp ( t) and using 
method of harmonic balance, one gets two sets of homogene- 
ous linear algebraic equations. As all the Fourier constants 
can not vanish simultaneously the coefficient matrix must 
bo singular. This condition loads to two equations between 
E and q in implicit form. From these one can finqdtransi- 
tion values of E and q using a Nag-subroutine C05AJF [17] 
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which attempts to locate the zeros of a continuous function 
by a continuation method using a secant iteration. For the 
function, given in the form of a determinant, obtained in 
the present case, the determinant is calculated by use of 
another Nag-subroutine F03AAF [18] which utilizes crouts 
factorisation method, 

3 • 1 ANALYSIS OF INSTABILITY REGIONS : 


A standard Mathieu equation which governs a wide 
variety of problems in mathematical physics is 


d^w 

dz 2 


+ ( § +E Cos z) w = 0 


For the system considered the corresponding equations for 
two mode approximation are 

EK f 2 


q B l,zz + E B l,z + ^ 4 
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40B 2 n □ n 
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- 4q M 22 ) Cos z B 2 


ol 


.. (3.2) 
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These are two coupled differential equations which 
are obtained after linearising the equations (2.23 and 2.24) 
obtained by applying Galerkin’s method to the governing 
partial differential equation of motion with two mode 
approximation. As the coefficients in these equations 
have periodicity of 2% one has 2n periodic and also 4n perio- 
dic solutions of these two equations for obtaining the 
transition points from stability to instability. 

(a) For 2n Periodic Solutions: 


For 2 7i - periodic solutions one can write the 
solutions of Bj. and Bg by the following Fourier series 



+ 


N 

2 

n=l 



+ 


N 

l 

n=l 


(a n Cos n z + b n Sin n z) 
(a n Cos n z + Sin n z) 


....(3.3) 


Substituting the equations (3.3) for B^(t) and (t) into 
equations (3.1) and (3.2) and comparing the harmonics, one 
gets recurrence relations for a n » b n ? a / ! and • In the 
present analysis, N is taken as 3. In fact if N is increased, 
one can get more accurate results but that is at the expense 
of higher computation time. Having done the harmonic balance 
one gets a system of simultaneous, homogeneous algebraic 
equations where the unknowns are the Fourier constants. Write 
this system of equations into a matrix form 


[CM2] 


14X14 


[ a 0 a i b i a 2 b 2 a 3 b 3 a 'o 


*1 b| a^ b£] 


" ^1X14 


(3.4) 
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whore CM2 is tho coefficiont matrix for 2n periodic solution 
given on page no 30 . Since all the Fourier coefficients 

can not vanish simultaneously for tho non-trivial solution, 
the coefficient matrix has to bo singular i.e. 

det [CM2] =0 .... (3.5) 

This is an explicit equation in E and q. To get 
the transition values of E and q set the value of E equal 
to a prescribed value and solve equation (3.5) for q. For 
E = 0, one gets the transition points on q axis, called as 
bifurcation points wherefrom the instability boundaries 
bifurcate. For this the solution procedure is as follows: 

Start with E = 0 and any one of the values of 
q obtained therein, increase E by some amount and look for 
solution of q from the matrix equation (3.5) with an initial 
guess of q using the Nag subroutine C05AJF [17 ]. ( It 

attempts to locate a zero of a continuous function by a 
continuation method by using a secant iteration). To locate 
the transition curves in tho neighbourhood of a specific 
bifurcation point q^ _q , search for values of q for a 
proscribed B by initial trial values of q^ - sA q and 
q E _0 + q. Continue increasing E and obtain the correspon- 
ding values of q in the neighbourhood of qg__Q by adopting the 
current value of q as the initial guess. Repeat this 

procedure for all the bifurcation points q (i) ,i=0, 1, 2 . . . . 

E =0 
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( b ) For 4 tc Per iodic Solutions : 


D 


= 2 (a n Cos §5. + b „ Sin nz 


n= l , 3 , . . . 


n 


N 


For 4% periodic solution of equations (3,1) and 
(3.2), one 'writes the following Fourier series solutions 

2 ^ 

... (3.6) 

3 0 = 2 (a' Cos —■ + b ' Sin ) 

n=l 3 ^ 

Substituting (3.6) with N = 5 in equations (3.l) 
and (3.2) and comparing the harmonics the matrix equation 
obtained is similar to that of 2%— periodic solution and 

is given by £ CM4 J 22 ^ 22 1- ^1 ^’3^*3 ^ 5^5 ^1^1 ^3 ^3 ^*5 ^5 

= 1X12 (3,?) 


where [CM4] is the coefficient matrix for 4n periodic solution 
given on page no 31 . This matrix is singular as equation 

(3.7) is a system of homogeneous simultaneous algebraic equa- 
tionsin Fourier coefficients which do not all vanish simul- 
taneously for non- trivial solution. So 


det [ CM4 ] = O 


(3.8) 


It is an explicit equation in E and q describing 
the transition values B and q for 4% periodic solution of 
(3.1) and (3.2). The method for obtaining the transition 
values has been discussed earlier. 
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3.2 RESULTS : 

Firstly, the principal parametric instability 
region is found out. To do this, one mode ( the first one) 
approximation of lateral vibration is made i.e. putting N=J 
in the first equation of (3.6) and neglecting the second 
equation of (3.6). Therefore, a set of homogeneous simulta- 
neous algebraic equations is a 1 and bj^ is obtained i.e. in 
equation (3.8) determinant [CM4] is of order 2X2 from which 
E versus q graph is plotted for the following values of the 

F o 

parameters d, = 0, M^,. = 0.1, 50 and p— = O . The 

cr 

results obtained ( Fig. 3a) are comparable with those obta- 
ined by Gurgoze [14], For the same one mode approximation 
and taking three terms of the Fourier series, the determinant 
[CM4] in equation (3.8) becomes 6X6. Again the principal 
parametric instability region is found out as shown in Fig. 3(a). 
Here a slight change in the instability boundaries is noticed. 

A' 

The bifurcation curves have also been obtained for 6^ =0,1 
by taking N=1 and N=5. It is observed that the instability 
boundary is lifted up for both the cases of N. This is 
similar to one reported by Gurgoze and Nayfeh. 

Results for two mode approximation have also been 
obtained and are shown in Fig. 3(b). The principal parametric 
instability region obtained is also compared with that of 
Gurgoze for one mode approximation. 
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Frequency ratio,q(= 2 ^) 

Fig. 3a Principal parametric instability region for one 
mode approximation (compared with that by 

G'urgoze) with M*=0.i, K, =50 and-^— =0.0 

r * fcr 
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CHAPTER 4 


METHOD OF HARMONIC. _BALA N.CE 


The method of harmonic balance is to express the 
periodic solutions of the differential equation (D.E.) 
of motion by Fourier series (F.S.) , substituting this in 
the D.E. and equating the coefficient of each of the har- 
monic to zero to get a system of algebraic equatiorsrela- 
ting the frequency and the amplitudes of different harmonics 
of the system response. The accuracy of the resulting 
solution depends on the value of the amplitude of first 
harmonic and the number of harmonics assumed in the solution. 

In this chaptey, the steady state response for one 
mode approximation at the principal parameteric resonance 
are obtained, one similar to that of Gurgoze [14], To do 
this the F.S. solution is put into the governing D.E. of 
motion for ono mode approximation and the coefficients of 
Sine and Cosine terms in the resulting equation are equated 
to zeros. This results in two simultaneous homogeneous alge- 
braic equations in two variables { coefficients of Sine and 
Cosine terms in the F.S.) because the governing D.E. is homo- 
geneous and it has cubic non-linearities. Those two algebraic 
equations are non-linear in nature but by introducing the amp- 
litude of steady state vibration, one can make it a system 
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of two linear homogeneous algebraic equations. For a 
non-trivial solution of tho transverse vibration of the 
beam, the determinant of the coefficient matrix must be 
zero. That gives a quadratic equation in the square of 
the amplitude of the steady state vibration as a function 
of system parameters M*, K^, G and frequency ratio q. 

4.1 ANALYSIS; 


The differential equation (2.18) , obtained from 
one mode approximation of the latral vibration of the beam 
is rewritten as 


fjft) + ™ B(t)+ (n + y Cos tot) B(t) + ( B, B, B ) =0 


(4.1) 


whero, t = 


oo 

“ol - E K f2 ’ 




2 

~W w o + k f 4 » 


(|J. (B,f B*) 


c ^3 3 tt 2 c 
5 B _ -3— - 


™.p .t, 

a B 


+ M 46 b b 2 + M 44 b 2 b' 


(4.2) 


and the overdots denote derivatives with respect to time t. 

In the non-linear system, one obtains periodic 
solution with frequency equal to that of the excitation -- os 
well asfb/n, n=2,3,..,, called subharmonic oscillations. 


When the excitation frequency happens to be equal to that 
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to be equal to that of the* subharmonic, one experiences 
the phenomenon of resonances, called subharmonic resonances. 
Similarly, superharmonic oscillations are experienced. When 
the frequency of excitation is close to one half of the 
frequency of free oscillation term and the system response 
is very largo, it is called the principal parametric reso- 
nance. 


To investigate the steady state response in the 
neighbourhood of principal parametric resonance, that is 
w 2 u) q ^ assume the 4 it - periodic solution of equation 
(4.1 ) as 

B(t) = a^Sin ( ?p-) + b-^Cos (yp>) .... (4.3) 


Substituting (4.3) into (4.1) and setting the 
coefficient of Sin ( and Cos (“-j|) separately equal to 
zero the following equations arc- obtained. 


(l-3n#+ ^q 2 ) a-2c> d qb + a ^[ ( ) a-;-q Orb - Kq 2 a]=0 

1 o 1 J- 


3& 


(l-H* “4 m 2 q 2 )b i +2§ d qa+ i[ ( j~) b-q D L a “Kq 2 b] = 0 

^ol 


(4.4) 


1 ~ - ' 4w ^; 1 ^ “ L 1 _i t 


whore, a^and b^comprise the amplitude and phase, 


a 


= 'CTTb 2 


amplitude 




of steady state vibration, 
K f* 




37 


M* 



9 


m 1 = I (B m 22 - |) = e ( u 2 /2) (m*+ §) - | ' 


m 2 = | (B M 22 + -) = B (ti 2 /2) (M*+ | ) + | , 

2 2 

r-. 2 M 3 n 

K = — Lie M* - 3 ~ - ^ ] 


and = 


3t: 2 < 

TT b d 


(4.5) 


Equation (4.4) has obvious solution a =b^=a = 0» For a non- 
trivial solution consider equations (4.4) as a system of 
linear homogeneous algebraic equationsin * and ' b^ invol- 
ving 'a 1 and set the determinant of its coefficient matrix 
equal to zero. This gives a quadratic equation in a having 
the solution 


a = ([ ( m, g * + 2 S d q 2 D L ) ± [ ( * + 26 d q 2 D L ) 2 

- (%*+46 2 q 2 ) (S* 2 +q 2 D l 2 )]^ 2 ]/( S^+q 2 q 2 ) ) 1/2 

.... (4 * 6 ) 

where 


S* — Kq — 8^. 


L =■ 


$ 


'* " : ”2 = [ I + T2' ] / (1 ~p 


2 


1 . K f* 


4w 0l 

m*= 1—2 IT* — q 
and 1 %*= (1-h*) (1-3^) 


), 


% 


cr 


0 f 

+4q 0 ( 1 — ) — m 0 (l-3p.^) ] 
4 

-16 m 2 m 2 q 


.... (4.7) 
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Tho equation (4.6) is an explicit equation for ’a’ 
as a function of q, the freqency ratio. It is called the 
response curve. 


The backbone curves, which indicate the relation 
between the bending eigen frequency to and the amplitude 
’a'of the non-linear beam system are obtained, equating 
terms to zero which correspond to damping and displacement 
excitation. To do this, replace ^ by to in equations (4.1) 
and (4,3) and proceed in the same way as that for dynamic 
response. This is equivalent to setting = 6 =0 

in the equation (4.6) and getting 


a 


= t(l-q 2 )/(Kq 2 -S.)] 1/2 


(4.8) 


where q 


ol 


ratio of eigen frequencies of the non- 


linear and linear systems respectively, 


4,2 RESULTS: 

Backbone curves obtained from equation (4.8) are 

graphically represented in Figs. 4(a) and 4(b) for 

F 

(i) 6 =0 , M*= 0.1, = 0.5, £^=0 ancl ae a 

cr 

parameter ranging from 0 to 50 and 

F 

(ii) B = 0, K f *=15, ^ =0.5, and M* as a parameter 

ranging from 0 to 0.4, respectively. 

From equation (4. 6) the response curves are 
plotted for the following cases 
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(i) E =0.05, M* = 0.1, =0 , p~ = 0.5 and stiffness 

r cr 

ratio as a parameter ranging from 0 to 50 , See Fig. 4 (c) 

F 

(ii) 6=0.05, = 2 , S', = 0, =0.5 and mass ratio I\l ^ 

^cr 

as a parameter ranging from 0 to 0.4, See Fig. 4(d), 

F 

(iii) M* = 0.1, K^^=15, S , = 0, ?£■ = 0.5 and excitation 

F cr 

amplitude E as a parameter ranging from 0.0 to 0.05, 

See Fig. 4(c) . 

Here both the stable and unstable branches of response curves 
are found out. If (+) ve sign is taken in place of (+) 
sign in the equation (A. 6) stable branches are obtained and 
for unstable branches (-)ve sign is taken. The stable and 
unstable branches are shown by firm and dotted curves 
respectively. 



Non - uimen 
amplitude 
response 



Fig. 4a. Backbone curves for various stiffness ratios 

F 0 o 

with M*= 0.1, 6 =0, — - = 0.5 and 0 .= 0 



Non-dimensional amplitude 
of response (a ) 



Non-dimensional amplitude of response la; 



Frequency ratio, q( s ^) 


Fig.4c Response curves for various stiffness 
ratios and 6 =0.05, M* = O.l, °-=0.E 

6 * • o.o 



Non-dimensional amplitude 
of response (a) 
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Stable branch 



Frequency ratio, q (= 35^) 

Fig. 4d Response curves for various mass ratios and 
6 = 0.05, K,*=2 ,-^j =0.5,6 d = 0 



Non-dimensional amplitude of response (a) 



Fig. 4e Response curves for various excitation 

amplitudes and M*= O.l, K flk = I 5,-^- =0.5 
6 d =0 
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C HAPTER 5 

ANALYSI S BY METHOD OF SLOWLY VARYING PARAMETERS 

In this chapter, the steady state response 
characteristics in the neighbourhood of main parameteric 
instability i.e. w is found by the method of slowly 

varying parameters* The backbone curves are found out for 

w -“or 

For linear undamped free systems the solution can 
be written as a Sine or Cosine function of time with constant 
amplitude, phase angle and circular frequency. However, for 
weakly non-linear system with small damping the solution for 
the displacement and velocity can be expressed like that of 
the linear case but with a slowly varying amplitude and phase. 
’Slowly' means that amplitude and phase remain nearly constant 
during the period of oscillation. This can be viewed as a 
transformation from one dependent variable to two dependent 
variables namely the amplitude and phase. So from this 
transformation, one gets two non-autonomous differential 
equations in two dependent variables, amplitude and phase. 
These non-autonomous differential equations are transformed to 
autonomous system by taking average values of the right hand 
side over the time period. It is presumed that the solution 
of those two autonomous equations represent approximately 



the solution of the non-autonomous equations thus leading 
to the steady state response of the system. The backbone 
curve is found out in a similar way. 
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5,1 ANA LYSIS s 

For one mode approximation of lateral vibration 

of the beam the following second order non-linear homogeneous 

DE was derived in Chapter 2 (equation 2,21) 

p 


B + £— 
, zz pio 


B 


m i “*G Kj ^-Fo v -i “ 

+ [ — +e( — - Mp 2 ) Cos ZJ B 


to 


to 


2 2 
* % 


+ — + K^) B ^ B B, 


3tT 

~sr 


c 

u oiium 

(UtO 


— Q»«~« 


+ 


M 44* 


B 


zz 


(5.1) 


For linearised version of equation (5.1) with 
no damping force and parametric excitation,, tne solution, in 
the neighbourhood of main parametric instability region 
(u)'iJ2a) ol ), can be written as 


B = a Sin ( § + ty) 

Where ’a 1 and 1 tjj 1 are constants. But for using the 
method of slowly varying parameters the solution of equation 
(5,1) should be written as 

B= a(z) Sin ( | (z)] ,....(5,2) 

where amplitude * a 1 and phase ’l[J r are time dependent. 

As stated earlier it is assumed that the velocity 
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B' = | a (z) Cos [ | + l|j(z) ] ....(5.3) 

where subscript dash means differentiation with respect to 
new time variable z. 

Dir f erentiating (5.2) with respect to z and 
comparing it with (5.3) one gets 

a* Sin Q + a l|J* Cos Q = 0 ....(5.4) 

where, Q = | + ijj (z). 

Differentiating (5.3) with respect to z 

■§"= Cos ( |f(jj) - |( | +l|j») Sin (| + iji) (5.5) 


Substituting (5.2), (5.3) and (5.5) in equation (5.1) one gets, 


- ~ (a'CosQ-allJ' Sin Q) = - | Sin Q + 0 P — a o o~ 

l+M 44 a 2 Sin 2 Q 


w ol -BK„ ICn 

+ a [— -! + G ( ) Cos z] 

(0“ 0) 


Si nQ 


2 . 2 
TU (0 ^ 

+ “2 (”6~ + K f4 ) a 

( jO 


p 9 

1+M 4 4 a Sin Q 


3 Sin°Q 

— — I n 11^ - 

l+M^ 4 a 2 Sin 2 Q 




VI. X SinQCos 2 Q 

+ a d TTaX ” 13 

4 l + ^44 a Sin Q 


3* 2 c 


_2 




o Sin QCosQ 

, i Vf-najv . m * i m — i i^l . « i^i i^ m A 

9 9 

l +M 44 a ~Sin Q 


(5. 6) 
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Solving a* and iji 1 from equations (5.4) and (5.6) 


» _ 


2 a Sin Q Cos Q + a 

1+M 44 a Sin Q 




K 


+ 2a [ ^ W— + G (-42 - m ) Cosz] Sin^,Cos.O 


0)- w 2 *22 


l+M 44 a 2 Sin 2 Q 


2 2 

_ . _ -2 , y \ 3 Sin Q Cos Q 

w l+M 44 a^Sin 2 Q 


+ K ( 


^46 3 S inQ Cos°Q 

2 a ” 1+^44 a 2 Sin 2 Q 


3ti 2 c _3 Sin 2 QCos 2 Q 


.(5.7) 


and 


^ ^ Sin 2 Q + 


o 


c Sin Q Cos 0 
^ l + i;,,a" 2 SirAT 


44 

. , r “4- E K f2 + p , K f2 
+ 2 [_ 4- + E (. 


to 

2 2 

2 , 71 % 
+ ? ( T - 




Wj 2 ) Cosz ] — 


2„ . 2, 


l+M^a Sin Q 


+ K f 4 ) a 


2 


Sin ‘Q 

”” o 0~" 

l+i'/l 4 4 a ‘"Sin‘"Q 


■4.3, 


+ M 46 2 ^S-CosV. _ 3rdL g 2 

a 1+M 44 a 2 Sin 2 Q 8 ^ 1+ii^a 2 Sin 2 Q 


(5.8) 


Up to this point, no more than a co-ordinate transformation has 
been carried out and the differential equations (5.7) and 
(5.8) aro still exactly equivalent to (5.1). If the non- 
linear t^rms arc small then a 1 and iji’ aro also small, that is , 

and 

the amplitude 1 a'^/phase r IjJ * change only slowly. r Slowly' 
here moans that the value of l{j in the argument (j| +l|j) as well 
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as the value of a(z) remain nearly constant during a time 
interval of duration T 0 = 2 %/ ( = 4m. Since it is genera- 

lly not possible to obtain an exact solution to (5,7) and 
(5-8) those equations are simplified by replacing the 
right hand sides by their temporal mean over the interval 
[t,t+T Q ] « In forming the mean value, a and l(j arc kept 
constant on the right hand side of (5.7) and (5.8). So 
equations (5.7) and (5,8) are replaced by 


•a 


i 


i_ 

4 it 


An 

| / Sin 2(~H-ljj)dz 


4~ 


IJXjoa 


/ 


4tc Cos 


2(z +W 


4 tc w* 

+ a / [“- — 

o 


44 

-e k 


-w— 

l+M /lzl a 2 Sin 2 (|+(lJ) 


■ -dz 


f 2 


0) 


2 


L 9 

+ B( ifi 


Sin2(|+!(J) 
M 22 ) Cosz] A 


+ 


+ 


0) 


M 44 . 3 r 4Tt Sin(~ +l|j) Cos 2 (|fl|j) 

)^\**'*** cl J MMVM# ir^Vn.pnLnrlnT' WW <M » **— 


2 


, y )a 3 J 4% Sln 3 ( § 4) Cos(|+l(l) 
f 4 


l+iV^ 4 a 2 Sin 2 (| +l[j) 


dz 


8 


2 ~ r y / w ^2 

L l+M^a 2 Sin 2 ( ^4-fjj) 


44 c 
dz 


4m 


44 1 
2 tz 


l+JVi 44 a Si n 


2<-,4 “ 2 (§ +ty) 


elz 


3m' 


8 (.to) 


Sin 2 ( f +TjO Cos 2 (f +</) 
1+M 44 a2 Sin 2 (§+l]j) 


.. ( 5 . 9 ) 
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and 


r 


4n 


[-5 7 "sin 2 (| +l|j)dz + - c - f 4% Sin (§H[J) C os(f+l|)) 


4fC L 2 


2 


pw 


1 +M^ 4 a ^ Sin 2 +t|l ) 


dz 


4tu K__ . 2/z 1 . 

+2 f [ — — w — — + 6 ( ~H -A’ ) Cosz] Sin ^2 + ^ 

0 w w ^ 2 _. 2 


1+1VL a 2 Sin 2 (§ +l|j) 


dz 


44 


2 7I~C0 


/ 

+ -75 t 


2 2 
ol 


(0 


8 


+ W ^ 


471 Sin 4 (f»+l|j) 


1+M 44 a2sin2 (| ♦W 


dz 


+ 


M 

**++. 

2 

.2 


4 6 


2 j 47T Sin z (| +l|j) Cos*(|+l|l) 


, 2 /z 


2 /z 


1+M 4 ,a z Sin 2 (§ +ljj) 


dz 


o 


44 « wj.ii \2 

4rc r.j ^3 /2, 

2 


§{£ £_ 2 / TC Sin d (§+l|j) Cos (§■ 4) 

8 pm J — — —— — — — 


1+M 44 a z Sin 2 (| +l{j) 


dz] 


respec tively. 


5.10) 


All the definite integrals in equations (5.9) and (5.10) are 
calculated subject to the condition 


1 + M 44 a > 0 


where is negative for mass ratio less than 
0,37 and it is the coefficient of non-linear inertia of 
the concentrated mass. Therefore for a high value of 'a' 
the method becomes invalid. But the results can be obtained 
satisfactorily for a definite range of amplitude. For steady 
state solutions of the autonomous equations (5.9) and (5.10) 
l|j' = a* = 0 



c ,:--i , , . * * ~ , , a . , 

1 1 , 


A104Jk.2.n. 
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giving 


( M 44 Q 2 +2— 2 ^l\ 4A q 2 +1 ) ( - 2Gq 2 M 09 )Sin 2(|j 


M/ 


[ 44 a 2 [l~^ 4 a 2 +l + — - (i^ 4 a 2 +2-2 1/ AL 4 a 2 +l)]qS' 

1 Tt.T r ‘ r r ( 


3'/i 


16M 


44 


d 


(5.11) 


and 


( 2 ^ Vl 44 ^ +1 ~ M 44 a 2 ~2)(^*~2E q 2 M 22 ) Cos2llJ 

^4 A ^ 0~“ r ° 9 1 

4 v M 44 a +1 [q^« 7-~ (1- f— -5 ) (1-24,,) 

/4 44 a M 4/ r a +1 


1 


M 


44 


[ 1- 


s A 


M fl 2 ( 1+ Sr*"? rr )] ■tr" 2 

M 44 a M 44 a +J - 


J ol 


JVU f 

35P* [ 1+ 

*14 


2(l-^M 44 a 2 +l ) 


44 
"44 


3 q 2 3 


(5.12) 


where 




he- = 


_c 

2 pw 


ol 


EK 


203 


f 2 

ol 


s 


2 2 
tc oC 
0 

8 


+ K. 


f4 


and 


q 


CO 


2w 


ol 


Al = M., a 2 +2-2 ^ a 2 +l. 


Introduce 
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bi = Y7— T“ 

M 44 a +1 , 


and 


Cl = (n # - 2Eq 2 M 22 ).Al, 

D1 = M 44 a 2 [l-Bl + — .2s!] q S 


D2 


2 4 
‘44 


M44 *lF 


d 


M/| A Ci 


/P" Bl. [q* 


M 44 a ' 


Bl ^ 


4 [1 ~ (1 ' ^ )] 


4&j 


M 


46 A1 2 -I 
2 ~ 2 ^ J * 


Mp 

r rr 


Therefore, the equations (5.11) and (5.12) can be written as 


Cl . Sin 2l|j = D1 (5.13) 

and - Cl . Cos 2(jj = D2 (5.14) 

Equations (5.13) and (5,14) are equivalent to 

Cl 2 =D1 2 +D2 2 (5.15) 


from which the steady state response characteristics can be 
obtained. 

The backbone curves are obtained by substituting 

by to in equation (5.1) , assuming B=a(z) Sin [z+l|j (z)J 

and procoding in the same way. Finally the equations obtained 

are exactly similar to (5.11) and (5.12) but here q ~r~— . ■ 

ol 
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5.2 RESULTS : 

Equation (5.15) is an implicit equation of 
amplitude ’a' and frequency ratio ! q’ and it is very diff- 
icult to express 'a' as a function of 'q 1 . So first find 
out the value of q for a=G and then find its value in the 
neighbourhood of a=0. The procedure for computing q for 
different values of 'a' in the neighbourhood of a=0 is the 
same ns one adopted earlier for finding the instability 
regions, section 3.1(a). 

Based on the computed results for E =0, 5^=0 and 
F 

J2. = o, the backbone curves have been plotted for the follow- 
F cr 

ing cases: 

(i) Mass ratio M*=0.1 and stiffness ratio K f * from 0.0 to 
50.0, Soe Fig 5(a). 

(ii) Stiffness ratio = 15.0 and mass ratio M* from 0.0 
to 0.4, soc Fig 5(b). 

C F o 

The response curves for o d =0 and p— - 0.5 and 

0 

other parameters mentioned below have been shown: 

(i) Excitation amplitude E = 0.05, mass ratio i% =0.1 and 
stiffness ratio K f * from 0.0 to 50.0, see Fig. 5(c). 

(ii) Excitation amplitude E = 0.05, stiffness ratio K f *=2.0 
and mass ratio M* from 0.0 to 0.4, see Fig. 5(d). 
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(iii) Mass ratio IV^-O.l, stiffness ratio =15.0 and 
excitation amplitude G from 0.0 to 0.05, see Fig. 5(e). 




Non-dimensional amplitude of res 



Frequency ratio , q b (= ) 

Fig.5b Backbone curves for various mass ratios 
with G =0, K f = 15,-4“- =0.5 and S d =0. 0 



Stable branch 
-Unstable branch 


Frequency ratio, q ( = 

Fig. 5c Response curves for various stiffness ratios 
with 6 = 0.05 , M* = 0 . 1 , -£r=0.5 and S d = 0 
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Fig, 5d Response curves for various mass ratios 
with 6=0.05, Kf* = 2,"~p 0.5, 8 d -0.0 
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CHAPTER 6 

MALYSL S BY METHOD OF M.JLT.TPT.F. SCALES 

In this chapter, the method of multiple scales 
has been adopted to obtain the steady state response of 
the one mode approximation of the system. 

The method of multiple scales is to consider the 
expansion representing the response to be a function of 
multiple independent variables or scales instead of a single 
variable. The number of independent time scales needed 
depends on the order to which the expansion is carried out. 
Here, the expansion of the dependent variable l? (response) 
is carried out only upto first order. The equation of motion 
for one mode approximation is written as a modified Mathieu 
equation having damping and non-linear terms which are assu- 
med to be of order E ( non-dimensional excitation amplitude 
being considered as small). In this equation, the assumed 
expansion of tho response is put and the coefficients of like 
powers of G are oquated to zeros leading to a sot of second 
order ordinary differential equations. For a bounded solu- 
tion of these equations, one has to avoid secular terms by 
making the coefficient of secular terms to be zero. This 
gives a differential equation (dE) in amplitude ’A' from 
which tho response of the system can be obtained using polar 
form for ’A". The real and imaginary parts of this resulting 
equation are separated into two non— autonomous D.^is which 
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are converted into autonomous system by appropriately 
eliminating the time ' t*. From the resulting equations the 
steady state response is found out. 

6.1 ANALYSIS : 

Equation (2.21) involving independent variable z can 
be rewritten as 


B(t) -(- B(t)+[w 2 .,+E (K f0 ~M 00 w 2 ) Cos wt] B(t) 


J ol 


f2 22 


e f(B, b,b") 


.... ( 6 . 1 ) 


where t= z/m , C = c 

6 ' 


• • * 


f(B,B, B ) 


e 


[§b : 


B +M 46 B -B 2 


+ M 44 B 2 B- ] + K f2 B 


.... (6.1a) 


and overdots denote differentiation with respect to t. 


In the method of multiple scales, one introduces 
the nt?w independent variables T Q and according to 


T n = G n t for n= 0,1 


dl 


n 


.n 


dt 


....( 6 . 2 ) 


So (tiyfollows that the derivatives with respect to 
time t become expansions in terms of the partial derivatives 
with respect to according to 
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d _ 
d t “ 


dT a di 
o d u 1 

dt dT o + dt~ 



D 0 + GD x 



9 

dt 



+2E D q D 1 + E 2 D 2 


(6.3) 


The solution of (6.1) can be represented by an 
expansion of the form 

H(t,e) » B 0 (T 0 ,T 1 ) + B % ( T 0 > T x ) (6.4) 

Whore , T Q ~ t and = Gt. 


Substituting (6.3) and (6.4) into (6.1) and equa- 
ting the coefficients of G° and G on both sides, the follwoing 
equations arc obtained 

D 2 B + u) 2 , B =0 (6.5) 

o o ol o ' 1 


D o' V “ol h =~ 2D o D l V (K f2 -M 22 “ 2 ) COS cot B 0 

“ ? D o V f < 5 0 - 1 W D o (6 ’ 6) 

Equation (6.5) is homogeneous linear DE with constant coeffi- 
cient and (6,6) is non-homogeneous DE. 

The solution of (6.5) is written as 

0 O = a(T x ) exp ( i« 0l T 0 ) + Xuq exp (-i0> ol T o > 

(6.7) 
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where , ' A 
c ornp 1 o x 
term i a 
solu tion 


15 an unknown complex function of q and A is its 
coniugate. It is seen that 3„ is real as second 
complex conjugal 

5 0 


to or the first. Substituting the 


oi into equation (6,6). 






( K jr 2~* 2 ^ "* ) [ 


2i t0 o iA’ (Tj_) exp ( im ol T o ) 

a (t i } 

~ m ' (“ iw ol T o +i wT 0 ) 


4- oxp (~iu) ol T 0 “imT o ) ) ] 

- fr ico oi A(T i } ox p ^ iw oi T o )+ cc 

+ ;t '[A(f 1 ) oxp (i w 0 j_T 0 ) + A (T^) exp (~ioj ol T o ) , 
i %l^ A(T l ) 0X P (i V T o) “ A (Tj,) oxp (-ico ol T 0 )] f 
~ (jJ ol ^ A ( Tj ) Gxp l*» 0l l 0 )+K h ) exp (-iw ol T o ))] (6.8) 


where CC denotes tho complex conjugate of the preceding terms 
and primes denote cli ff orentiation with respect to Ti. 

. To investigate tho steady state response in the 
neighbourhood of principal paramo tor ic resonance write 


oo = 2 oj q1 + Bcr 

whore a ~ detuning parameter. 


( 6 . 9 ) 


^(w-Woi) T 0 = w ol T o +aT l (6.10) 

Mow from the equation (6.8) the solution of can 
be looked for. Tho particular solution of this equation has 
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some secular terms, containing the factor T exp (+ ico T ) 
and in order to get a bounded solution of B x and, therefore, 
a uniformly valid expansion of the response, these secular 
terms arc to be suppressed. Making the coefficient of this 
secuJ ar term as zero 

A 


A ' 


£ i %1 A + 


(l\f 2 - ^ exp(icrT 1 ) 

2ti/ w 0 j_ 

*~r ( f exp(-iu ol T Q )dT^ 


2 it o 

O 


.... ( 6 . 11 ) 


whom f(M 0 ,VV D o ®o^ ec I ua ^i° n (6.8) is expanded in a 
Fourier nurles [7] according to 


o 


f <V D o'\,’ D o U o> = n l 0 ^Wo-W exp(in %1 T 0 ) + CC] 

and thus the coefficient of exp (iw ol T 0 ) in it is 

2%M q1 


CO 


fx(U 0 - D 0 15 0 - U P 0 )= # / f CS 0 ,D 0 B 0 , Do 5 o> ex P (-i“olT 0 ) dT 0 


As 'A' is a complex quantity it can be expressed by 

.... ( 6 . 12 ) 


I 

Ap- a exp (ire) 


where n and a aro real. 

Putting the equation (6.12) into (6,11) and separating 

the real and imaginary parts, the lollowing two equations in 

'a* and ’«* aro obtained 
„ , a (K. FO ~M 09 o) 2 ) Sin (aT,-2a) 


4w 


2ji 


ol 


3 


r 27C / W Ol 

L- J f . Sin(w , T +a) dT 
2 tc o ol o c 
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and 

act' - (K f2 - ^22 w ^ ) Co s ( cj T^-2a ) 

2%/b) , ( 6 < 13 ) 

•MW I V -X- 

2tc £ f .Cos (w ol T 0 +a)dT 0 

Putting (6.12) into (6.7) , one gets 
:3 0 » a Cos (^ 0l T 0 +a) 

~ a Cos 0 .... (6.14) 

Where 0 « ai ol T o +Q: 

Therefore, the zeroth order approximation for B is 

given by 

B = a Cos 0+0 (G) ....(6.14a) 

The equations (6,13) are non-au tonomous as time T^ 
appears explicitly in these equations. Substituting 

l|j - oT^ - 2a 

in equations (6.13) and thus converting these into an autono- 
mous system of D.E.s 

a ' ~ ^ K f2~ M 22 w ^ Sin ^ ** 2p 

, 27t 

- — / f. Sin 0 60 

2mo ol 0 
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oil 1 ao a ,, o. . 

2 = 2“ ~ (K f2 -M 22U ^)Cos l|j 


1 2111 

2tcw , ^ f. Cos 0 d0 

01 o 


where (jj* a 
2 ~ 9 


2 ~ a ♦ 


. . . (6.15) 


For its steady state solution 
a* =4' s=0 


. . . (6.16) 


Put tho zeroth order approximation of 3 (given by 
equation (6,14a) into (6.1a), integrate the definite integrals 
in equations (6.15) and use the condition (6.16) to get 


(K f2 _M 22" 2)Sin ’I' = ?k 


(6.17) 


3s~ (K f2 -,Vl 22“ Z) Cos * = syy [ 


3o M 46 _3M 


46"”° 44 ,2 

“TT - -“olJ 


Eliminating Ip from equations (6.17), one gets an equation 
for amplitude 'a' of the response as a function of the detuning 
parameter c and the amplitude of the excitation £, called the 
f requncy response equation. On multiplying the resulting 

I-*.- Zr.2 /, ,2 N — f* t /-l ir* o 


equation by (G 2 /w^) on both sides 
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( EK f2 


2\2 / c s2 


4u> 


r = (■, 


Ol 

BK 


■■ 22 - ' - 


+[ — !*-■ + q-l 

n. 


2to 


ol 


where, 


q= " 2 “ol 

c _ 3 S' 

~ “o' 

4 “ol 


( 5 * - K )] 2 


1 + , 


2m 


ol 


and K 


3Al 44~ M 46 


For zero damping (c=0) 
a = [ 


. M*+q-i + ( ?r - e M 2 2 q2 ^ 1//2 


| (S*-K) 




i. . . (6*18) 


(6.19) 


The backbone curves are obtained from equation (6.1) 
by replacing ^ by oo , following the above procedure with 
w=w o i+Ea and equating the terms to z'-ro which correspond to 
damping and displacement excitation. This is equivalent to 
setting c=E=0 in equation (6.18), where q is to be interpre- 
ted as : 

w ol 


E 1 


_q-l_ 


(S*-K) 


] 


1/2 


( 6 . 20 ) 


The expressions for d* and K are same as those in 
the equation (6.18). 
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6.2 R ES ULTS: 


The backbone curves are plotted for the following 
two cases from the equation (6.20) 

F o 

(i) B =0, M* = 0.1, p-— » = 0.5, = 0 and as a parameter 

ranging from 0.0 to 50.0 ( See Fig. 6a) 

F 

(ii) 6=0, K^. = 15, = 0.5, 6^=0 and i\. as a parameter 

r cr 

ranging from 0.0 to 0.4 (See Fig. 6b). 

Tho response curves are plotted from equation (6.19) 
for the following cases taking only (-) sign in place of (+) 
sign in it 

F 

(i) G=0, 05, M# = 0.1., o, =0, = 0.5 and stiffness 

h cr 

ratio Kc. x . as a parameter ranging from 0.0 to 50.0 ( See Fig. 6c) 

(ii) G = 0.05, K-* = 2 , £ ,=0, ^ = 0.5 and mass ratio as 

^cr 

a parameter ranging from 0.0 to 0.4 (See Fig. 6d) 

F 

(lii) M* = 0.1, K f * = 15, 6 d = 0, p — =0.5 and excitation 

0 X 

amplitude E as a parameter ranging from 0,00 to 0.05 ( See 
Fig. 6o). 



Non-dimensional amplitude of response ( a ) 



Fig. 6a 


Backbone curve for various stiffness 


ratio with mass ratio M* = 0. 1, € = 0 
^ - 0 . 5 , 6,-0 

Vl 



Non-dimensional amplitude of response (a) 


i 



Frequency r 
Fig.6b Backbone curve for vari< 

Ku- 15, e -0,f; -0.! 


0.4 


0.2 

0 .! 

M*= 


0 



Non-dimensional amplitude of response (a) 



Frequency ratio, q (=- 55 ^) 

Fig. 6 d. Response curves for various mass ratios 
with 6 =0.05, K f ^= 2, p^"=0.5, 6^=0 
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CHAPTER 7 


DISCUSSIO NS 


The principal parametric instability region 

for the idealised engine-valve mechanism for mass ratio 

= 0.1 f stiffness ratio = 50 y ratio of the preload 

in the spring to the beam critical load «-£ = 0.0 and 

F cr . 

for one mode approximation with no damping (o ) is shown in 
Fig* 3a. For small values of S , the transition curves are 
linear and comparable with that by Gurgoze[i4j but for higher 
values of Q there is a side tilt from the linear transition 
curve. For the same system parameters and one mode approxima- 
tion with damping 5^ =0*1? it is observed that this region is 
lifted up as was noticed by Gurgoze[l4]. Considering three 
terms in the Fourier series repres antation of the response 
and taking the same system parameters f this region is relo- 
cated for <*^=0 .0 and 0*1 . Comparing the transition curves 
for one term and three terms of the Fourier series expansion 
for B( t), it is observed that the right side droops down with 
marginal change in the left side. Taking two mode approxima- 
tion the transition curves starting from q = 0.5 , UO, 2.0 
and other admissible discrete values equal to and less than 
4 are plotted in Figc(b)with the same system parameters and no 
damping. It is observed that the principal parametric region 
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is reduced and the other instability regions are narrow . 

The backbone curves found by the method of 

slowly varying parameters (SVP) are plotted in Fig. 5(a) for 

constant values of the system parameters = 0.1 , G = 0, 

F 0 

P^" r ~ ? a d an< 3 stiffness ratio varying from 0.0 

to 50.0. Similar curves are found for the system parameters 

kf# ^ = ®> —2~ = 0.5, 5^ = 0 and various values of M # 

(0*0 ~ 0.4), It is observed that the nature of the change 
in these curves for different system parameters by this 
method is similar to that found by application of harmonic 
balance (H.B) method by Gurgoze and also by the application 
of method of multiple scale (M.S.). The response curves 
(stable and unstable branches )f ound out by method of S.V.P. 

F n 

for the system parameters B = 0.05, = 0..5, -0, 

c r 

M# = 0.1 and various values of ( 0.0-50.0) are shown in 

Fig , 5(c) . Similarly response curves are shown in Figs. 5(b) and 
5(e) to show the effect of system parameters M* and S . It 
is observed that the nature of these curves is similar to 
that found by the method of M.S. in Chapter -6 and also to 
that found by the method of H.B. done by Gurgoze . Backbone 
curves and response curves found by these three methods are 
compared in Figs.7(a ~e). 



Norrdimensionai amplitude of response(a) 



Frequency ratio, q b = (zr^) 

Fig. 7.a Comparison of Backbone curves obtainec 
by three different methods with M|O.I, 

6 = 0.0, 6 d =0.0 and -^ = 0.5 
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j. 7b Comparison of Backbone curves obtained 
by three different methods for Kf*= 15.0 

6 = 0.0, 6(j= 0.0 and = 0.5 
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Frequency ratio, q (= 2S57 ) 

Fig. 7c Comparison of response curves by three different 
methods for M*= 0. 1, € = 0.05 , 6 d = 0.0 and -Jr- =0 
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CHAP :ER 8 

CONCLUSION 

The non-linear parametric vibrations of a 
restrained vertical beam with a concentrated end mass at 
the top under harmonic ground excitation has been studied. 

The instability region(s) of the linearised system for one 
mode and two mode approximation is/ are found out. The 
steady state response curves are obtained by the method of 

(i) harmonic balance, 

(ii) slowly varying parameters 
and (iii) multiple scales. 

The major findings are summarized below. 

The principal instability boundaries bifurcating 

from — ~ — = 1 are not linear for higher values of the 

2w ol 

excitation amplitude. The effect of the second mode on 
these boundaries is observed to reduce the instability region. 

For a fixed excitation frequency 

(i) an increase in the stiffness ratio decreases 
the amplitude of the steady state response 

(ii) an increase in the mass ratio increases the 
steady state amplitude and 

(iii) an increase in the amplitude of the ground 
excitation increases the steady state amplitude. 
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APPENDIX -I 


Combining equations (2,1) and (2,2) 

- F(t > v-Mj+Hj-M, a EXv ,,(l~v 2 7 1/S (1.1) 

Expand the right hand side of the above equation and keep 
only upto third order principal non-linear terms to get 

- F(t) v-M 1+ M 2 -M 3 =EI V (SS (1+ | v^ s ). (1.2) 

Differentiating the above equation twice with 
respect to s 


/ d 2 M 3 d 2 M „ d 2 M„ 

— Fit) v — — -.T . ~ ^ ~ 

,ss ds 2 ds 2 

~ ^- V ,SSSS + 2 V ,SSSS V ,S +3V ,S V ,SS V ,SSS + V ,SS ] • 

(1.3) 

From the equation (2.9) 

M 1 = / Hv, tt [ /(l- 5 v 2 ) d ] dp 

S S 7(1 

L 

- R / (1- | v 2 p ) dp (1.4) 

In equation (2.7) M-q is defined as 
L P i ,2 

M 11 = / P v tt ^ ^ 1 ”2‘ ,T] ^ dT ^ d|3 ( x * 5 ) 

s 1 s 

L 

= / f(P»s)dp, 

s 

P 12 

whore f(p,s) = pv ^ / (1- g v ,g ^ dri 

^ s 
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Differentiating partially with respect to s by Leibnitz 


rule 


c) " 

- / f(p,s)dp 
s 

L v 

~ / a7" f (i 3 > s ) dp - f (s,s) 

s 

L p 

s ds"~ ^ v ,tt ^ 2 v Jr] ^ dr| ^ ^ 

s 

- >‘ v ,tt { S d- k\l ) dn 

L a 

= { ^ v ,tt [ { “ (1 " I v ,l > d i 

- a- 1 v^s) + a_| v 2 ? ) a| ] d p 


{ d - 1 v i’s ) d P 

s 

( 1_ 2 v ,s ^ ^ L ^ lv ,tt 
s 


( 1 . 6 ) 


Again differentiating — ’ — - partially with respect to s 

c)s ’ 


( 1 - | ) 


ST { t* v ,tt d P 


+ v ,s v ,ss { d P 


^ v ,tt (1 - k v ,s ) + v ,s v ,ss { ( 1 ' 7) 


5 \s' T \s \ss % i",tt 
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Therefore, one gets 


_ d \i a 2 ro , l i 2 

as 2 " as 2 "" ~ Ts”" 2 ' [R { {1 ~ 2 v J dp 


“Att (1 - 2 v f s ) + v ,s v ,ss A d P - R \s \ 


( 1 . 8 ) 


[As ’R 1 is constant ] 


From the equation (2.10) 

„ r l r J , o P 

fvU - J p fg + + + f (v . + v „ v , . ) dn !.(/ a 


,nt + v ,n v , n tt ) d, i 3-d v^dn) dp 


Differentiating M n partially with respect to s 


(1.9) 


dM 2 r L P 

t*“— ei: / 8 r / /»r 


{ Is ^ l(g -tt + / V .« + +) d 8) 


»T)t s T) j T) 1 1 ' 


dn) ] dp 


{ (v fnt +v ,v V ,vtt } aR] d'v^dr,) 

y O 


L °) 

=- V jS / [H(B tt + /(V * t -w V tt ) dn] dp (1.10) 

s o 

aw 2 

Again dif forontiating partially with respect to s 


jpr “ ^ v ,s [g ,tt + { (v ,nt +v ,r) v ,ntt )dT)J 


L Jo 

r r , r / £ 


■^.SS { [9,tt + { ( V ;„ t +V (i V )tltt ) dn] dp 


(1.11) 
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From the equation (2*11 ) 

1 .2 


L J 


ig ~ I cv >t U U- | )cJn] dp 


. ... (1.12) 

s » u s ' 

Observe that can be obtained by substituting c for p 
and v .j. for v ^ in the expression for (given by 
equation (1.5) ). Therefore, substitute c for p, v ^ for v ^ 
and Mo for A/L^ in equation (1.7) to get 


d 2 M, 

V 

— 

ds -1 ' 


cv 




(1- k v 2 J + V V 


ss 


L 

/ CV 


1 1 


dp 


(1.13) 


Substitute equations (1.8), (1.11) and (1.13) in 
the equation (1.3) and simplify to get the equation (2.14). 
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